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ABSTRACT 

The numerical solution of heat and mass transfer thermal diffusion effect on MHD free convection flow of 
stratified viscous fluid. Here we considered Visco-elastic, Darcy resistance terms, the constant permeability of the medium 
and neglected induced magnetic field in comparison to applied magnetic field. The velocity, temperature and concentration 
profiles are presented graphically and we obsereved that the temperature does not fluctuate with change in), 
K (Permeability parameter), Gr (Grash of Number), A (Thermal diffusion parameter), and Sc (Schmidth number). 
The concentration profile increases in increase in A But decreases with increase in Sc. Velocity profile rasises with 
increase in Gr, k, and A but it decreases with increase in M(Magnatic field). The equations are solved by Galerkin finite 
element method. 
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INTRODUCTION 

Heat transfer is that science with seeks to predict the energy transfer which may take place between material 
bodies as a result of a temperature difference. Thermodynamics teaches that this energy transfer is defined as heat. 
The science of heat transfer seeks not merely to explain how heat energy may be transferred, but also to predict the rate at 
which the exchange takes place under certain specified conditions. The fact a heat transfer rate is the desired objective of 
an analysis points out the difference between heat transfer and thermodynamics. Mass transfer is the net moment of mass 
from one location to another, mass transfer occurs in many process i.e., such as absorption, evaporation, drying, 
membrance filtration and distillation. Mass transfer is used by different scientific disciplines for different process and 
mechanisms. The phrase is commonly used in engineering for physical process that involves diffusive and convective 
transport of chemical species with in physical systems. In heat and mass transfer the transfer of heat is between a surface 
and moving fluid (liquid or gas), when they are at different temperature which is known as convection. 

In geothermal reservoirs and geothermal we can find very important application of convection problem in a 
porous medium. And also we can observe the process of heat and mass transfer has application in chemical process 
industries, fluid fuel nuclear reactor, aeronautics and many engineering applications in which the fluid is working as 
medium. The large verity of technical and industrial applications has stimulated enthusiasm on the study of heat and mass 
transfer in the convection flows. Although the effect of stratification of the medium on the heat elimination process in a 
fluid is important, very little work has been recorded. 
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Ostrach[1953] studied free convective flow past a vertical plate. Even siegel [1958] investigated transient free 
convection from a vertical flat plate. Cheng and Lau [1977] and Cheng and Teckchandani [1977] obtained numerical 
solutions for the convective flow in a porous medium bounded by two isothermal parallel plates. However, porosity 
measurements by Benenati and Broselow's [1962] noticed that porosity is not constant varies from the surface of the plate 
to its interior, to which as a result permabulity also varies. Soundalgekar [1972] with regard to unsteady free convective 
flow stuided the effects of viscous dissipation on the flow past an infinite vertical porous plate. Chen et al [1980] have 
studied the combined effect of buoyancy forces from thermal and mass diffusion on forced convection. Bejan a and khair 
[1985] have investigated the mass transfer to natural convection boundary layer flow driven by heat transfer. The free 
convection on a horizontal plate in a saturated porous medium with prescribed heat transfer coefficient was studied by 
Ramanaiah and Malarvizhi [1991]. Lin h t and wu c m [1995] considered the problem of combined heat and mass transfer 
by laminar natural convection from a vertical plate. The mass transfer effects on free convection flow of an incompressible 
viscous dissipative fluid was studied by Manohard and Nagarajan [2001. Rushi Kumar b and Nagarajan [2007] have 
investigated mass transfer effects of MHD free convection flow of an incompressible viscous dissipative fluid past an 
infinite vertical plate. Recently, the effect of stratified viscous fluid on MHD free convection flow with heat and mass 
transfer past a vertical porous plate was discussed by Agrawal [2012]. 

In the present section we have considered the problem of Agrawal by introducing thermal diffusion and solved by 
Finite element method. 

MATHEMATICAL ANALYSIS 

The study of Heat and Mass Transfer Thermal Diffusion Effect on MHD Free Convection Flow of Stratified 
Viscous Fluid with assumptions of temperature of plate is constant, Viscous and Darcy's resistance terms are taken into 
account with constant permeability of the medium, The suction velocity normal to the plate is constant 

(V' — —U 0 ), Boussinesq's approximation is valid. 

We considered a system of rectangular co-ordinates 0 (x ,y ,z ) where as y — 0 on the plate and leading edge 
is z axis. The density variation with temperature is considered while of all fluid properties remain constant. Due to the 
influence of density variation in the terms of the momentum and the energy equation and the variation of the expansion 
coefficient with temperature is neglected. The variations of density (p 0 ), viscosity (pi 0 ), elasticity (ct 0 ) and thermal 
conductivity (k 0 ) are considered as 

p = p 0 e- b 'y', n = Mo e- b 'y', a = a 0 e -» V, k T = k 0 e -» V 

where p 0 jU 0 , cr 0 , and k 0 are the coefficients of density, viscosity, elasticity, and thermal conductivity 
respectively at y = 0, b' > 0 represents the stratification factor. 

Under these conditions, the problem is governed by the following system of Equations: 

Equation of continuity: = 0 (1) 
Equation of Momentum: 

p S + v t) = U"%) - K + £ ) ul + - T '-> + ° (2) 
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Equation of Energy: ^- + V' ^- — — 4~ (k T (3) 

M 6J df dy' pC p dy'\ ' dy'J v ' 

Equation of Concentration: + V = D 7-77 + Di ^-yr (4) 

M df dy dy 2 1 dy 2 ' 

Where, u , v are the velocity components, T , C are the temperature and concentration components, v is the 

kinematic viscosity, p is the density, o is the electric conductivity, B 0 is the magnetic induction, k T is the thermal 

conductivity and D is the concentration diffusivity, Cp is the specific heat at constant pressure, Dl is the thermal 
diffusivity. 

The boundary conditions for the velocity, temperature and concentration fields are: 



u = 0,T = T w , C = C w oty = 0 
u =0,7 = T ml C = C^aty' -» 00 

Introducing the non-dimensional variables 



v J u 0 (Tw-rj (c 1 w -c 1 0O )' 



(5) 



t = — -,y — ^—^-,K — ^~2~ is the permeability parameter, U 0 is the reference velocity, 



M 



= is Magnetic field parameter, P r — is Prandtl number, 



Sc = - is Schmidth number, G r — vg ^ Tw Tc °- ) j s Modified Grashof number for heat transfer, 

D r Ul 

v4 = 1 w " is the thermal diffusion parameter , N 0 = — - — is the buoyancy ratio, 

(C w — C (») p (7w — Toe) 

/? is the thermal expansion coefficient, p 1 is the concentration expansion coefficient and, b is the stratification 
parameter. Other physical variables have their usual meaning. 

Introducing the non-dimensional quantities describes above, the governing equations reduce to 

0 (6) 



9V 
~dy~ 



( l - fo )S = 0 _ ( M+ 9 u+G '- (0+iVoC) (7) 



du 

dt v ~ J dy dy 
dd . b . 90 _ 1 d 2 9 

Tt~( 1 ~V r ) d-y'-V r ^ W 
dC dC _ 1 9 2 C 9 2 0 

^ _ 97 _ s79y 7 9y 2 ^ 
and the corresponding boundary conditions are: 

(10) 

y — 7 j 

METHOD OF SOLUTION 

By applying the Galerkin finite element method for equation (7) over a typical two-nodded linear 



u = 0,0 = 1,C = laty = 0 
u — 0, 9 = 0, C = 0asy->co 
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element (e) ( y ; - < y < y k ) is 

£[£-^(0-»>^-^-(«+3^*»-)]*-° 

where R* = G r {9 + N C C),N = [Nj,N k ],0^ = 
U W = w . 0Wf N . = ^f,7V fc = ^J, IW = y k - Vj = h, 
The element equation given by 

y ; k#; kJ k^/ jv k jvj kJ y kty w fc jvj kJ ay + 

/ i\\NjNj NjN k -\rUj-\ \Njl 

Now put row corresponding to the node i to zero, from equation (13) the difference schemes with — h is 
^ Wl _ 1+ 4 il + ^ + ((M + i)^«-i)u I . 1 + e + (M + i)au I + (( M + i)^«4)u wl = 



(11) 

(12) 



(13) 



R and here R = R 



(e) 



(14) 



Using the Cranck-Nicolson method to the equation (14), we obtain: 

AiU 1 ^ + A 2 ul +1 + A z u[H + A A u\_ x + A s u{ + A 6 u J i+1 = R (15) 

Bjg + B 2 e> +1 + b 3 8^ + B 4 eU + b 5 o( + e 6 e/ +1 = o (16) 

dc/L^ 1 + C 2 c/ +1 + C 3 c/^ + Qi t + C 5 c/ + Qc^ = R 1 (17) 
^ = (l - 6r + 3p(l - b) + (m + fc) , i4 2 = ^4 + 12r + 4 (m + ^ fc) ,,4 3 = (l - 6r - 3p(l - b) + (m + ^ kj, 

A 4 = (l + 6r - 3p(l - b) - (m + kj , A 5 = (4 - 12r - 4 (m + kj ,A 6 = (l + 6r + 3p(l - b) + (m + kj 
B 1 = (l - 6r £ + 3p(l - bj) , B 2 = (4 + 12r ,B 3 = (l - 6r £ - 3p(l - ft)), 
£? 4 = (l + 6r^-3p(l-6)),B 5 = (4-12r^), B 6 = (l + 6r^-+ 3p(l - 6)), 
d = (l-6r^+3p), C 2 = (4 + 12r^),C 3 = (l-6r^-3p), 
C 4 = (l + 6r^-3p),C 5 = (4-12rj-),C 6 = (l + 6r^-+3p), 



Ri = ( 



0i +1 -2^9i+9i_ 1 
h 2 



) 
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Here where h=0.1 and k=0.05 mesh sizes along y direction and time direction considered for computations. 

A t Xi = B u i = 1,2,3 .... (18) 

Where Ai's are matrices of order n and X t and B t are column matrices having n-components. 

RESULTS AND DISCUSSIONS 

Velocity distribution of fluid flow is plotted and shown in figure 1 for n=0.1, Pr=0.71, No=1.5, t=0.1, b=0.1 for 
with different value of M, Gr,A ,K, and Sc. figure 1 shows that all velocity graphs are increased sharply for y=1.2 as such 
velocity in each graph begins to decrease and tends to zero with the increasing in y. It is also observed from figure 1 the 
velocity decreases with increase in M but it increases with the increase in Gr, A, K and Sc. But the temperature remain 
constant parameters taken for velocity. The concentration profiles for different values of diffusion parameter A are shown 
in figure 2. It is noticed that concentration profiles increases with the increase of A the concentration profiles are higher 
when Sc is 0.4, compared with Sc is 2. It is also observed that concentration increases with the increase in A, but it 
decrease with the increase in Sc in figure 2. The graph of Sc verses the C presented in figure 3. The Sc effect is leads to 
decreases the concentration profiles are observed form the figure. 

CONCLUSIONS 

• The velocity increases with the increase in A (Thermal diffusion parameter). 

• The concentration also increases with the increase in A. 



i 

0.5 



Figure 1: Velocity Profiles 




Sc=2.0. A=0.2. 0.4. 
Sc 4:' A n: 0.4. 0.6 



Figure 2: Concentration Profiles 
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Figure 3: Concentration Profiles for Different Values of Sc 
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